The phase shift effective range expansion from supersymmetric quantum mechanics 
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Supersymmetric or Darboux transformations are used to construct local phase equivalent deep and 
shallow potentials for t partial waves. We associate the value of the orbital angular momentum 
with the asymptotic form of the potential at infinity which allows us to introduce adequate long- 
distance transformations. The approach is shown to be effective in getting the correct phase shift 
effective range expansion. Applications are considered for the ^Pi and partial waves of the 
neutron-proton scattering. 
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I. INTRODUCTION 

In the description of the interaction between composite 
particles by local potentials an ambiguity arises between 
different phase equivalent or isophase potential families. 
There are shallow potentials, which possess only physi- 
cal bound states and deep potentials, which possess both 
physical and unphysical bound states. The latter, called 
Pauli forbidden states, simulate nonlocal aspects of the 
potential, or else, the complexity of the interaction be- 
tween composite particles. The number of Pauli forbid- 
den states can be predicted from a microscopic descrip- 
tion of the interacting particles The application of 
the supersymmetric (SUSY) quantum mechanics f5| to 
the inverse scattering problem provides an elegant and 
powerful algebraic way to understand the relation be- 
tween such phase equivalent potentials 0, 0, IE 0- A 
supersymmetric transformation can be seen as a specific 
Darboux transformation. In the following we shall use 
Darboux and SUSY transformations as synonyms. 

The above procedures cannot directly provide a correct 
behavior of the phase shift at energies small relative to 
the potential strength, i. e. a correct effective range 
expansion for £ > 0. In our opinion the reason is that 
the role of the angular momentum for a given central 
potential was not properly understood so far. In the 
framework of SUSY quantum mechanics the problem has 
been raised by Sukumar 3] and Baye and Sparenberg Q 
but not solved in principle. However, it has been tackled 
pragmatically in Ref. Q. For example, in the case of 
the £ — 4: partial wave it led to an S-matrix containing 
powers of k restricted to n > 10. How one arrives at 
such a restriction is neither explained nor the power 10 
is justified. 

The basic idea of this study is to associate the angular 
momentum with the long-distance asymptotic behavior 
of the potential, irrespective of its singularity at the ori- 
gin. This is in the spirit of Ref. 8\ where these asymp- 
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totic limits are independent of each other. This starting 
point will provide a new possibility for getting a correct 
effective range expansion of the phase shift which is the 
following Taylor series expansion in the vicinity of k = 
(see e. g. 0) 



r,2<' + l 



1 



1 



cot Se{k) = h -rof/e + 



dot 



(1) 



Here age is the scattering length and roe the effective 
range. The expression implies that for a given £, in 
the series expansion of tan Si{k) the coefficients of the 
terms containing powers of k below 2£ -\- 1 must vanish. 
In the frame of SUSY quantum mechanics we solve this 
problem by introducing adequate long-distance Darboux 
transformations. 

The paper is organized as follows. In the subsection 
A of the next section we introduce the Darboux trans- 
formation method and briefly review the ^-fixed transfor- 
mations. In subsection B we introduce the long-distance 
transformations. Section HI is devoted to results and ap- 
plications to the neutron-proton scattering. Details are 
worked out for the £=\ and £=2 partial waves. Conclu- 
sions are drawn in the last section. 



II. THEORY 

A. ^-fixed Darboux transformations 

We recall that the Darboux transformation method 
consists in getting solutions ip of one Schrodinger equa- 
tion 



hiip^Eif, hi = -— + Vi{x) 



when solutions ip of another equation 



hoii = El/;, ho 



Voix) 



(2) 



(3) 



are known. This is achieved by acting on ip with a dif- 
ferential operator L of the form 



-d/dx + w{x) , 



(4) 
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where the real function w{x), called superpotential, is de- 
fined as the logarithmic derivative of a known solution of 
(0 denoted by u in the following. One has 



w = u'{x)/u{x) , Hqu ~ au , 



(5) 



with a < Eq, where Eq is the ground state energy of ho 
if it has a discrete spectrum or the lower bound of the 
continuous spectrum otherwise. The function u is called 
transformation or factorization function and a its factor- 
ization constant or factorization energy. The potential Vi 
is defined in terms of the superpotential w as 



Viix) = Vo{x) -2w'{x) . 



(6) 



Eq. Q defines a first order Darboux transformation. In 
the following we shall deal with chains of N successive 
transformations of this type. 

Let us start by first considering £- fixed transformations 
as in Ref. P|. This means that we use a special chain of 
N — 2n~v first order Darboux transformations with v > 
0, generated by the following system of transformation 
functions 



Vl{x) ,Vy (X) , Uy+l {X) , Vy + l (x) , 

hoUj{x) — —a'juj{x) , hoVj{x) 



,Unix) ,Vn{x){7) 

-b]v,{x), (8) 



where Vj are regular [vj (0) = 0) and Uj irregular {uj (0) ^ 
0), the latter being expressed in terms of the Jost solu- 
tions as 



Uj{x) = Ajf{x, —iaj) + Bjf{x,iaj) 



(9) 



They have arbitrary eigenvalues 



and —hj respec- 



tively, but always below Eq. If we are interested in the 
final action of the chain only, the solution 'il)]\r{x, k) of the 
transformed equation with the Hamiltonian 

hN = S/dx^ + Vn (10) 

corresponding to the energy E ~ is given by [Tol | 

ipN{x,k) = W{ui, . . . ,UN,i>o{x, k)) W^'^{ui, . . . ,un) 

(11) 

where W are Wronskians expressed in terms of Uj, de- 
noting symbolically any function of Q and of '>pQ{x,k) 
which is a solution of the original Schrodinger equation 
corresponding to the same energy E. In the Hamiltonian 
()10|l the transformed potential is 



,Un) 



(12) 



For iV = 1 one has W{ui) = ui and one recovers © 
with u = ui. If Vq is finite at the origin, Vjv behaves 
as I'lv + l)x~'^ when j; — > 0. Therefore the parameter 
h' is called the singularity strength. The formulas 
and H12|l result from the replacement of a chain of N 
first order transformations by a single TVth order trans- 
formation, which happens to be more efficient in practical 
calculations. 



In Ref. la] we obtained that the transformed Jost func- 
tion Fjv is related to the initial Jost function Fq by 



FNik)=Foik) Yl 



i = i 



k + ibj 



n 



ibi 



(13) 



For h' — the first product is unity. Since a Jost function 
is analytic in the upper half of the complex A;-plane (see 
6- g- UM) 9-11 b's must be positive whereas the a's can 
have any sign, so that every positive Oj corresponds to a 
discrete level E — —aj of /iat. 

The corresponding phase shift (k) can be written as 



Sf{k) = S^{k) + Af{k) 



(14) 



where 5f{k) is the initial phase shift due to the potential 
Vq and A^(fc) is the phase shift produced by the chain 
of N Darboux transformations 

n n 

^?'(^) ~ ^ arctan(fc/aj) — arctan(fc/6j) . 

J=U + 1 J = 1 

(15) 

This is consistent with the asymptotic form of the scat- 
tering solution sm{kx — + 6f). In the limit A; — > oo 
one has 6f ^ {£ — i^)7r/2, in agreement with Ref. for 
a singular potential of parameter ly. More detailed dis- 
cussion ofproperties of £-fixed transformations may be 
found in Q. 

In the case £ = by expanding S^k) and the arc tan- 
gent functions in IjlSfl in power series one obtains the 
effective range expansion For £ > the situation is 
more subtle, since the first term in power series of arc 
tangent functions is proportional to k. Therefore in the 
usual practice based on the SUSY approach, where S° {k) 
is fixed, it would be difficult to cancel the undesired pow- 
ers of k in order to comply with As mentioned above, 
we believe that the reason is that one deals with Darboux 
transformations which do not affect the long distance be- 
havior of the resulting potential, as it was pointed out in 
Ref. 13. 



B. Long-distance Darboux transformations 

To change the long distance behavior of a potential 
by SUSY transformations we use transformation func- 
tions with zero eigenvalue. To cancel undesired powers 
in the series expansion of arc tangent functions we derive 
a proper background phase shift as shown below. 

Consider the potential Vq (x) which for a; — > oo behaves 

as 

Vo{x) = ^^^t +0(x-^-^) , 7>0, ^>0. (16) 

As it is known jlll | the Schrodinger equation containing 
a potential satisfying (|16|l has zero eigenvalue solutions 
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with the following asymptotic behavior at x — > c» 

vix) = Cx^+^ [l + 0{x-^)] (17) 
uix) = ^[l + Oix--')] . (18) 

The functions (|17|l are regular at the origin but singular 
at infinity and the functions H18|) are just the other way 
round. When these functions are taken as transformation 
functions the change in the potential for sufficiently large 
X has one of the following forms 

AV{x) = -2[lni;(.T)]" = ^^^i^ + o(a;-2-7) (19) 

9f 

AV{x) = -2[lnM(a;)]" = - + 0(^-2-7) (20) 

where AV{x) ^ Vn - Vq with = 1 in It is clear 

from here that the function (|17|l increases the value of 
£ by one unit and the function (|18|l decreases it by one 
unit. Moreover a linear combination of (|17|) and H18|) 
is a function of type H17|) . They form a one-parameter 
family, while the function H18|) is uniquely defined (up 
to an inessential constant factor). In this family there 
is only one function regular at the origin. This func- 
tion, used as transformation function in the Darboux al- 
gorithm, changes both v and i but all the other mem- 
bers of the singular family change £ without affecting v. 
In the following we shall use the singular functions of 
the one-parameter family defined above to derive phase 
shifts leading to a correct effective range expansion. We 
shall therefore show that the parameters appearing in the 
linear combination of (|17|l and (|18|l can be chosen such 
as the resulting phase shift provides the general effective 
range expansion (Q. Hence, starting with a given Vq with 
£ = we first perform a number N = £ of transformations 
which give the correct long distance behavior of the po- 
tential Vg and introduce £ parameters in the phase shift. 
Next, an ^-fixed chain is performed, producing the final 
phase shift H14|) for which the potential Ve plays the role 
of the initial potential. The latter transformation does 
not affect the asymptotic form of the potential Vg at large 
X. Hence, the resulting potential Vn has an asymptotic 
behavior corresponding to the £th partial wave. The ad- 
dition of £ zero-energy eigenfunctions to the N = 2n — v 
first order transformation functions used in the ^-fixed 
chain increases N hy £ units. This means that in the A^th 
order transformation to be used below the total number 
of transformation functions is 



N ^2n 



£ 



(21) 



If we start with the zero initial potential, Vq = 0, the 
formula (|14|) for the phase shift has to be modified as 
follows: 



(22) 



Here (5^(fc) is produced by the long-distance transforma- 
tions which give rise to an intermediate potential Vg. In 



the following (5|(fc) will play the role of a background 
phase shift. The additional phase shift A^~'^(fc) corre- 
sponding to the ^-fixed subchain oi N — £ = 2n — v trans- 
formations has the same form as (|15|l . We shall illustrate 
this procedure by applications given in the following sec- 
tion. 



III. APPLICATIONS 

A. The case £ = 1 

We start with the potential Vq = in Eq. (0). Let us 
take zio,i = x + xq as the transformation function with 
zero eigenvalue, where xq > is a free parameter. Then 
the first order transformation operator Q takes the form 



dx x + xq 
The transformed potential is 

2 



{X + XqY 



xo > 



(23) 



(24) 



and its Jost solution fi{x, k) may be found by applying 
the operator H23|l on the Jost solution fo{x, k) = exp{ikx) 
of the free particle equation. After dividing by the factor 
ik one finds 



/i(x,fc) = 1- 



1 



ik{x + Xq) 



cxp{ikx) 



(25) 



For xo ^ the potential H24(l has = 1^1 = and its 
Jost function Fi{k) coincides with fi{Q,k). If one now 
applies the operator H23|) on an oscillating solution of the 
free particle equation sm(kx + Si) one obtains 



ipi{x, k) = ^kcos{kx + Sl) + 



sin(fca; + Sl) 



X + Xq 

This solution is regular at the origin provided 
6l = arctan kxo , 



(26) 



(27) 



and has the asymptotic behavior sm{kx + Sl ~ ^) at 
X ^ 00, it describes the £ = 1 scattering state of the 
potential (|24() . Thus we have switched from the partial 
wave f = of Vo = to £ = 1 of the potential Now 
we can perform ^-fixed transformations. The free param- 
eter Xq will be chosen so that the final phase shift 1)22(1 
will have the correct effective range expansion. Replac- 
ing S^ by the result H27|) and expanding all arc tangent 
functions in power series one can see that the coefficient 
of the term proportional to k in tan S^ vanishes for 



(28) 



3 =1^+1 
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Both the regular and irregular solutions corresponding to 
the potential H24|l can be found with the help of the Jost 
solution H25|l . But in the spirit of Ref. ^3 , we can avoid 
this step, thus considerably reducing the amount of nu- 
merical work. This means that in the formulas Ijlll) and 
()12|l we can directly use appropriate solutions of the free 
particle equation which are simple linear combinations of 
exponentials 

ipj{x, k) — Aj eKp(ikjx) + Bj exp{—ikjx) , Re kj = , 

where we have to find the correct ratio Bj/Aj. Since 
the regular solutions of the potential (|24|) . consisting of 
functions Vj{x) = Liipj(x, —ibj) satisfy the condition 
Vj{Q) — 0, which fixes the ratio Bj/Aj, we have free 
particle solutions of the form 

ipj{x, —ibj) = {bjXo + 1) exp{bjx) + (bjXa — 1) exp(— b^x) . 

(29) 

The irregular solutions for the same potential, defined as 
Uj = Liijjjlx, —ittj), should be obtained from the func- 
tions 

ipj{x, —iaj) = Aj exp(aja;) + Bj exp{—ajx) , (30) 

which for Bj/Aj ^ {ajXo — 1)/ (ajXo + I), Aj ^ and 
Uj > have an increasing asymptotic behavior but if 
Aj — they decrease asymptotically. To stress the dif- 
ference between solutions with different asymptotic be- 
havior we choose in the latter case aj < and Aj ^ 
but Bj = (for more details see 0). Here aj, bj and 
Xq are parameters of the model to be found below. Then 
in the A^th order transformation we use the functions 
uo,i = x + xq, and to calculate l(T^ . Note nev- 
ertheless that in the resulting phase shift, given by (|22|l . 
Sg has to be replaced by Si of H27|) since the initial poten- 
tial for the subchain of £-fixed transformations is now Vi 
of (|24|) . For the iVth order transformation the potential 
Vi and the phase shift 6l play an auxiliary role. 

As an application we look for neutron-proton (np) po- 
tentials which reproduce the 'pruned' phase shift of [l^l 
for the ^Pi partial wave. This phase shift together with 
the theoretical values obtained from the expression H14() 
and denoted by 6l are exhibited in Tabled The six fitted 
S'-matrix poles are 

ai = -0.7290, 02 = -0.7295, a^ = 1.0368 . . 
bi = 0.4403, 62 = 0.4408, 63 = 3.3818 

in fm^^ units. The superscript 7 carried by the phase 
shift is consistent with the formulas H21I) and H22|l and 
it implies that we used a 7th order transformation ac- 
cording to ifTTl) and (fT^ . Then from Eq. one gets 
Xq — 3.0578 fm. Now we can expand all seven arc tan- 
gent functions appearing in (|14|l in power series. This 
leads to a correct effective range expansion given by 

k^ cot (5[(fc) = -0.3182 - 3.1511fc2 + . . . (32) 

from which one can extract the scattering length aoi and 
the effective range tqi defined according to In Ta- 
ble ITU these values are compared with another theoretical 



TABLE L The ^Pi phase shift. The theoretical value is 
calculated from Eq. 1141 with N = 7 transformation functions, 
Eiab ~ "^"^"^p Ecm, Ecm ~ where ji is the reduced mass, 
nip = 938!27 MeV and m„ = §39.36 MeV. The experimental 
values are from Ref. [l^ 



Eiai{MeV) 








5l{deg) 


14 




-4.1944 




-4.27887 


42 




-9.01021 




-8.91287 


70 




-11.99126 




-11.9844 


98 




-14.42546 




-14.5092 


126 




-16.60093 




-16.7024 


154 




-18.59407 




-18.6533 


182 




-20.42528 




on ^^ Q7 


210 




-22.09942 




-22.0187 


238 




-23.61771 




-23.4936 


266 




-24.9813 




-24.8577 


294 




-26.19215 




-26.1263 


322 




-27.25323 




-27.3113 


350 




-28.16843 




-28.4228 


TABLE IL Theoretical values of the scattering 


; length and the 


effective range. 








Q.01 








Reference 


3.143 


- 6.302 


- 2.224 


21.976 


present work 


3.023 


- 6.895 






M 


2.736 


- 6.449 


- 1.377 


15.027 


Reid93 ([15] 


2.4 ± 1.3 


-12.6 ± 2.2 






m) 



model |]J. They are surprisingly close to each other. 
The phenomenological Reid93 potential 0| also gives 
similar values Moreover the scattering length is lo- 
cated in the interval deduced from a partial wave analysis 

In order to construct potentials giving rise to the phase 
shift 5\, the poles H31(l have to be associated with trans- 
formation functions defined by (|29|l and H30|l . The poles 
bj correspond to regular solutions of Vi {x) resulting from 
l|29|l . The poles ai and 02, which are negative, correspond 
to decreasing functions of the form H30|l with Bj = 
(j=I,2). It remains the pole 03 which is positive. If we 
take iJs/As ^ (03X0 - I)/(a3a;o -I- I) in Eq. we 
obtain from (|12|l a one-parameter (1^3/^3) family of one- 
level isophase deep potentials with the discrete level at 
E— —a\. But if we choose _B3/yl3 = {azx^ — l) / {a^XQ + 1) 
the initially irregular function moves into the regular 
family, the level disappears and one gets the uniquely 
defined shallow representative of the family of isophase 
potentials, denoted by V7. Fig. ^ shows this potential 
from which the centrifugal term has been subtracted. 
The potential V7, is quite close to the Reid68 poten- 
tial (13, represented in the same figure. Figure |21 shows 
one of the deep isophase potentials corresponding to 



[ A 

0.2 0.4 0.6 0.8 a;(fm) 

FIG. 1: ^Pi np shallow potentials. W (solid line) together 
with the Reid68 potential |]3 (dashed line). 
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S -0.1 
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-0.3 




0.5 1 1.5 x(fm) 



_4 L . — . . — , — . — . . — , — . — . . — , — . — . . — lJ 

2 4 6 X{fm) 

FIG. 3: Asymptotic behavior in natural logarithmic scale of 
the value of the ^Pi np potentials. The full line represents our 
shallow potential of Fig. 1, the dashed line our deep potential 
of Fig. 2 and the dotted line the updated flst phenomenolog- 
ical Reid93 'T^l potential, including one-pion exchange. 

back to the nucleon-nucleon interaction in each case. The 
fact that the shallow and deep potentials are asymptot- 
ically the same is physically correct. As mentioned in 
the introduction, deep potentials reflect the composite- 
ness (quark structure) of the interacting particles in the 
overlap region but, once the particles are well separated, 
they could be treated as point particles as in OPE the- 
ories, which means that beyond some distance the deep 
and shallow potentials should coincide. 



FIG. 2: ^Pi np deep potentials. The full line corresponds to 
the potential which gives the correct effective range expansion 
1321 1 and the dots represent the Reid93 potential ,13,|. 



— 0.351. This constant has been adjusted to get a 
potential as close as possible to Reid93 in the inter- 
val 0.4 fm < a; < 2 fm. This deep potential possesses a 
Pauli forbidden state of energy E = — a| — —44.46 MeV. 
Contrary to the Reid93 potential which is deep but finite, 
our potential behaves at origin as — 2/a;^. If necessary, it 
can be regularized as for example in Ref. . 

It would be interesting to analyze the asymptotic be- 
havior of our calculated potential to see if it is compatible 
with modern phcnomcnological potentials constructed in 
the spirit of meson theory, i.e. which include one-pion- 
exchange (OPE) contributions. This is precisely the case 
of the Reid93 potential, which includes OPE with neutral 
and charged pions (for details see Ref. Ts'l). The com- 
parison given in Figure |31 shows that the shallow poten- 
tial V7 and its deep partner are practically identical in the 
asymptotic region and that they are both extremely close 
to the asymptotic form of the updated Reid93 potential 
Thus our potentials are in the excellent agreement 
with expectations from Yukawa's OPE theory. We re- 
mind that we subtracted the centrifugal barrier from our 
potentials in order to make this comparison feasible. In 
solving the Schrodinger equation this should be added 



B. The case i = 2. 

Now we have to apply two subsequent transforma- 
tions associated with functions corresponding to E = 
and then, as above, a subchain of ^-fixed transforma- 
tions. After the first transformation with the function 
uo,i ^ x + xq, the potential Vi of has u — l/(a;-f'a::o) 
and u — {x + xo)'^ as linearly independent solutions at 
E — 0. Their linear combination uq 2 = cu + u, which is 
the transformation function for the second transforma- 
tion step defined by the operator 



-d/dx + u'q2{x)/uq,2{x) 



(33) 



contains two free parameters xq and c. These can be 
chosen such as the series expansion for tan(52(fc) starts 
at . The intermediate (or background) potential 



V2 = -2[ln uo,2(a;)uo,i(2:)]" 



Q{x + xo)[{x + xaf -2c] 



[c+{x + xqYY 

(34) 

obtained from the L = L2L1 Darboux transformation 
(for more details see 0) plays now the role of the initial 
potential for an .^-fixed subchain of transformations. The 
background phase shift (modulo tt) corresponding to V2 
is 



5r, — arctan 



3x0 — fc2(xQ + c) 



(35) 
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Note that the function ip2{x,k) = £2^1 sin(fca; + 5|) is 
regular at the origin and describes an unnormahzcd scat- 
tering state for £ — 2. In the formula (|22|l we have to 
identify 5^{k) with of (|35l) . After expanding in 

power series all arc tangent functions we find that the 
coefficient of the term linear in k vanishes for a;o given 
by ^ and 

n 

^=HE(V' + V')] (36) 

ensures the cancellation of the coefficient of k^ in the se- 
ries expansion of tan (5^. Now, to find solutions of the 
free particle equations, giving rise to the regular fam- 
ily of the potential (|34|l . to be used in (|12|) . we need 
eigenfunctions of — —dP/dx^ satisfying the condition 
L2Liipj{xj —ibj) = at a; = 0. They are given by the 
following linear combination of exponentials 

tpj{x,—ibj) = [3xo + 3bjxl + b'^{c + xl)]exp{bjx) 

— [3x0 — SfcjXg + 5^(c -f Xq) ] exp{—bjx) . 

(37) 

The irregular family still results from H3U|) subject to the 
condition that the ratio Bj/Aj is different from that pre- 
sented in l|T7|l . 

With a fit of a similar quality to that performed for 
^Pi we could reproduce the ^I?2 partial wave phase shift 
of with the following four poles of the S'-matrix 

ai = -0.2047, a2 = -1.9800 , 
bi = 1.2305, 62 = 4.9631 

in fm^^ units. From Eqs. I|28|) and H36|) we get x^ ~ 
-4.375 fm, c 116.08 fm^. This leads to the following 
effective range expansion 

k^ cot 6l{k) = 0.4496 -I- 10.9878/0^ + . . . . (39) 



where the superscript N = 6 represents 4 transformation 
functions associated with the poles H38|) plus two zero- 
eigenvalue functions mo,i and uo,2 defined above. This 
is consistent with the formula (|21|) with n = 2, 1/ — 
and £ — 2. From the expansion H39I) one can extract the 
scattering length 092 and the effective range ro2 defined 
in Eq. (1). These values are shown in Table ITU They are 
comparable to those obtained for the potential Reid93 




IV. CONCLUSIONS 

By working out these two particular cases we have 
shown that a new insight emerges into the role of the 
angular momentum of a central potential. If associated 
with the long-distance behavior of the potential, it al- 
lows us to introduce transformations that bring i free 
parameters in the background phase shift (Eqs H27() and 
(|35|l 'l. When a correct effective range expansion is re- 
quired, each extra unit of angular momentum imposes a 
new constraint on the whole system of parameters of the 
model such as the number of constraints coincides with 
the number of parameters in the background phase shift. 
For the particular cases oi £ ~ 1 and £ — 2 explicit so- 
lutions of the constraint equations are given. Thus, the 
extension of the method to ^ > 2 is straightforward. 

Before ending we should mention that the generalized 
Levinson theorem Q is always satisfied in our approach. 
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